Dimensionless cosmology 



Ali Narimani * Adam Moss ^and Douglas Scott ^ 

Department of Physics & Astronomy 
University of British Columbia, Vancouver, BC, V6T IZl Canada 



Abstract 

Although it is well known that any consideration of the variations of fundamen- 
tal constants should be restricted to their dimensionless combinations, the literature 
on variations of the gravitational constant G is entirely dimensionful. To illustrate 
applications of this to cosmology, we explicitly give a dimensionless version of the 
parameters of the standard cosmological model, and describe the physics of both Big 
Bang Neucleosynthesis and recombination in a dimensionless manner. Rigorously de- 
termining how to talk about the model in a way which avoids physical dimensions is a 
requirement for proceeding with a calculation to constrain time- varying fundamental 
constants. The issue that appears to have been missed in many studies is that in 
cosmology the strength of gravity is bound up in the cosmological equations, and the 
epoch at which we live is a crucial part of the model. We argue that it is useful to 
consider the hypothetical situation of communicating with another civilization (with 
entirely different units), comparing only dimensionless constants, in order to decide if 
we live in a Universe governed by precisely the same physical laws. In this thought 
experiment, we would also have to compare epochs, which can be defined by giv- 
ing the value of any one of the evolving cosmological parameters. By setting things 
up carefully in this way one can avoid inconsistent results when considering variable 
constants, caused by effectively fixing more than one parameter today. We show ex- 
amples of this effect by considering microwave background anisotropics, being careful 
to maintain dimensionlessness throughout. We present Fisher matrix calculations to 
estimate how well the fine structure constants for electromagnetism and gravity can 
be determined with future microwave background experiments. We highlight how one 
can be misled by simply adding G to the usual cosmological parameter set. 

*anariman@phas.ubc.ca 
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1 Introduction 



The study of the Universe at early times and on large scales might allow us to discover 
physics beyond the currently accepted models, namely ACDM for cosmology and the 
standard model of particle physics. Since cosmology probes energies considerably 
beyond those attainable by man-made accelerators, there is hope that we can find 
evidence for new physics in this arena. One speculative idea for extending our view of 
physics is to consider the possibility that the fundamental constants are not actually 
constant. Variation of fundamental constants, if realised, would play a key role in 
theories beyond the standard model, in particular for ideas inspired by string theory 
or extra dimensions. In such theories the coupling constants can appear as fields which 
evolve in time (see e.g. [2 [2]). This idea is most familiar today through the many 
studies inspired by claims of a redshift dependence in the fine structure constant (e.g. 

miiEiE]). 

Since the pioneering articles of Dirac in the 1930s [3 18] regarding variation of the 
gravitational constant, G, there has been a large body of work on the possible variation 
of several different constants (see e.g. [9]). Over time physicists have become more 
aware of the fact that only variation of dimensionless quantities can be meaningfully 
measured and discussed [101 1111 1121 113) . This is because any measurement in physics 
can be reduced to dimensionless ratios of quantities - in other words, a measurement 
of something is always made relative to some other quantity of the same dimensions. 
In fact this underlies the utility of "dimensional analysis" [141 1151 116[ as a physics 
tool. 

It is reasonable to ask if one can express the results of any general physical obser- 
vation in terms of a few dimensionless parameters. In most of cosmology the relevant 
physical quantities for specific problems are formed from the set P = {c, h, eo,e, G, rrip, 
mc,fc,T}, plus extensive variables, such as distance, mass, rate, redshift, etc. Then, 
from dimensional analysis, the only dimensionless quantities which can be constructed 
are 

Gml me „ kT 
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and combinations of them. Here Og is the "fine structure" constant for gravity (a 
notation which appears to have originated with Silk [T^), T is the background radi- 
ation temperature, and the other quantities are familiar physical constants [18]. The 
importance of the dynamical variable T, temperature, in this set, as opposed to the 
other dimensional constants, will be highlighted in the following section. 

This analysis shows that if a specific quantity is dimensionless (which should be the 
case if one is talking about physical measurements), and contains one of the members 
of P, then the other elements have to appear in a way that leads to at least one of 
the dimensionless ratios introduced in ([T}. 

There are many papers in the literature (e.g. [Tgi[^[^[^[^[M[^[^[77l[251 
1291 1301 13l1 1321 133[ ) which try to answer a question such as: "What would happen to 
observable X, if the speed of light or the gravitational constant or Planck's constant 
had a different value". Such questions are not well defined, since one can tune the 
other parameters of the dimensional set P, such that the dimensionless ratios in ([1]) 
remain the same. Since any observable is essentially dimensionless, it can only depend 
on these dimensionless ratios, and one cannot discuss what happens to the observable 
X if, for example, G changes. On the other hand one can meaningfully consider how 
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X might change if Og varies. In this paper we will explore consequences of adopting 
this dimensionless thinking to measuring observables in physical cosmology^ 

As a prelude to the rest of the paper, let us consider Big bang Nucleosynthesis 
(BBN) and the final abundance ratio of helium atoms relative to baryons - a mea- 
surement which is clearly dimensionless. We will discuss BBN in a dimensionless 
manner in more detail later on, but it may be helpful first to explain the logic and 
sketch the main idea. When we say "nothing happens to observable X if G changes" , 
one might naively think that changing G would, for example, change the expansion 
rate of the Universe and hence lead to observable consequences. The point is that 
there is no observable which depends solely on the expansion rate. There also has 
to be some other rate (like the recombination rate, neutron decay rate, etc.) that 
one is considering in the problem. The ratio of these two rates (the expansion rate 
and the other relevant rate) is dimensionless, and therefore has to depend on the 
dimensionless ratios and not just G. Based on this, the general dependence of the 
primordial helium abundance Yp — Yp (c, h, eo, e, G, mp, rric) can always be reduced to 
Yp = Yp (og Therefore, the dimensional constants should not be regarded as 

independent, and any variation in these constants should always be understood as a 
variation in the dimensionless ratios. 

In this paper we will consider different aspects of cosmology in relation to funda- 
mental constants. The next section deals with some difficulties which arise when one 
tries to track a variation in the fundamental constants within the standard framework 
of cosmology. In sections 3 and 4 we try to express the results of BBN and recombina- 
tion in terms of the dimensionless ratios defined in equation Q. We keep everything 
explicitly dimensionless in these sections, but also simplify some of the equations. 
They therefore lack the full accuracy required for comparing with data, but neverthe- 
less demonstrate the main physics. In section 5 we work through the publicly available 
codes for computing cosmic microwave background (CMB) anisotropics, make them 
manifestly dimensionless, and present the results of a variation of Og or Qcm on the 
CMB power spectra. It is shown in this section that variation in the gravitational con- 
stant brings a complication because of the need to define a cosmological epoch. We 
also perform a Fisher matrix analysis to show more explicitly how one could obtain 
incorrect results by being careless in adding constants in a dimensionful manner. In 
section 6 we briefly discuss how this work could be extended to other areas of physical 
cosmology. 

In cosmology, Newton's constant, G, is particularly important. This is because it 
enters into the dynamics of the background, the growth of perturbations and large- 
scale geometric effects. There have been many studies of how one might constrain 
G/G using cosmological or astrophysical data (e.g. [35] and references therein), as we 
will discuss further below. It seems odd that studies of Qem and /i are usually careful to 
point out the importance of only considering dimensionless constants, while papers on 
G (or even sections in review articles) tend to ignore this entirely. We already discussed 
the importance of this issue in an earlier paper [36] , and we found that in general most 
of the published studies are unchanged if one replaces G with Qg. However, as we will 
see, there appear to be exceptions in some cosmological applications. It is for this 

*In a somewhat different context the importance of using dimensionless numbers in cosmology was 
stressed earlier by Wesson 1341 . 

tAnd also and ctw, as we will see later. 
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reason that we would recommend keeping things dimensionless as much as possible. 

At this point in a cosmology paper it is conventional to mention ones choice of 
units and parameters. We should point out that do not set c = 1 or G = 1, or 
make any similar selection. We also refrain from using any particular cosmological 
parameters, dimensional or otherwise. 

2 Fundamental constants in cosmology 

In cosmology one requires an additional set of parameters (as well as the dimensionless 
physical constants) to specify the model. This is a little different from setting the 
values of (dimensionless) physical constants, since the cosmological parameters are 
(at least statistically speaking) chosen from among a set of possible universes, all with 
the same physics and physical constants. Some parameters (or ratios) may be fully 
deterministic, but others could have stochastic values. Someday we might have a 
theory which tells us exactly the value for some of the parameters, but in the current 
state of cosmology, we consider them as free parameters, which are unknown a priori. 
The usual set consists of (at least): {pB, Pm, Pr , PA, Ho,To, A,n, . . .}. These are the 
energy densities in various components (baryons, matter, radiation and vacuum), the 
expansion rate, CMB temperature, and amplitude and slope of the initial conditions. 
Usually the densities are expressed in terms of the critical density (which involves G) 
to give the set {0.b, ^m, ^n, J^a}. Parameters such as the amplitude A and tilt n of 
the scalar power spectrum probably depend on fundamental constants, but in a way 
which is as yet unknown. 

The obvious question that arises here is "how do we define the cosmological model 
in a dimensionless way?" Some parameters are simple, for example the power spec- 
trum descriptors: we can convert from A to the dimensionless power at some fiducial 
scale (e.g. A^(A;o) used by WMAP [37], o"8 derived from galaxy clustering or Martin 
Rees' [M]); and the slope n is already dimensionless and so offers no difficulty. 
The real problems come from the fact that the other main cosmological parameters 
are epoch-dependent quantities. 

We can imagine a gedanken conversation with an alien civilization to help sort 
this out. To communicate information on the background cosmology, we would need 
to discuss what we mean by "baryons", "photons", "neutrinos", etc., and then give 
their densities in some way. However, we still have to deal with the fact that many 
of the usual parameters (the fis, H , T, etc.) depend on time; this is the issue which 
distinguished the 9 parameter from the others in Hence, even if we have estab- 
lished that we live in the same Friedmann model as our alien friends, we still have to 
determine whether we are observing that model at the same epoch or not. 

The best way to do this would be to agree on a fiducial period in the evolution 
of the Universe and then discuss where we are relative to that. An obvious choice 
is the epoch of matter-radiation equality (discussed in a related context in [39]), but 
there are plenty of other possibilities: when the Hubble rate is equal to a particular 
reaction rate; when the Thomson optical depth is unity; when matter has the same 
energy density as the vacuum, etc. Assuming that the Universe has fiat geometry (a 
clearly dimensionless statement), then we can give the values of the fis at the agreed 
fiducial epoch, together with one number to fix that epoch, say the value of 9 at 
equality. Then we only need to give the value of one parameter today in order to fix 
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the epoch at which we hve (making sure this is a dimensionless parameter of course). 
This could be any one of the fis today, or the value of Zcq, or the value of Hoto, or 60 
(= kTo/m-p(?). Anything which is changing essentially monotonically in time, and is 
dimensionless, will do (so Hoto is fine now, but useless billions of years ago, when it 
was hardly changing). 

This gives rise to two complications when describing the observable Universe, which 
are not there when one is discussing models of laboratory physics. Firstly, if one is not 
careful, then it is possible to effectively fix more than one of the parameters today, 
leading to inconsistent results. The second part is that different choices of the "what- 
time-is-it?" parameter are not entirely equivalent, because some contain a dependency 
on other physics parameters. We will show effects of this in section 5. 

Take for example the choice of either the CMB temperature To, or Hubble constant. 
Ho (made suitably dimensionless, using the Compton coupling time, or some other 
timescale involving atomic physics). These two are related to each other via Or and 
Og using the Friedmann equation 

° 45 ^nh^ ml ' ^ ' 

Therefore, either Ho or To is enough to determine today's epoch, and it is not con- 
sistent to use the Friedmann equation and treat both of these as free parameters. 
However, it turns out that it matters whether we choose To or Ho, since the pair 
(To, Ho) in our Universe would not satisfy the above equation in a universe with a 
different a,. We return to this issue in section 5. 



3 Dimensionless BBN 

Big Bang nucleosynthesis (BBN) is an area of astrophysics which has been thoroughly 
investigated for signs of variation in the fundamental constants (see e.g. [T] 1401 141j 
for non-gravitational couplings and BBN). Among the constants which are effective 
during BBN, the gravitational constant plays a key role, and hence there have been 
many published studies of the effects of a varying G on the primordial abundance 
of the light elements g2]|43|4l|43|4i[471[48l|49l[50l|5l]. However, since G is a 
dimensional constant, we calculate the abundance of helium synthesized during BBN 
in an explicit dimensionless way in this section, focusing on the dominant parts of the 
physics and neglecting some of the finer details. Our calculation, though crude, shows 
the role of Og in primordial nucleosynthesis and explicitly reveals the dimensionality 
in the relevant physics. 

The key parameter in BBN is the ratio of number densities of neutrons to protons, 
which can be defined as: 

— =e-". (3) 

The quantity u is the ratio of mass difference to freeze-out temperature: 

(mn-mp)c^ . 
kT, ■ 

Tf is explicitly the temperature at which the following reactions freeze out: 

e" -I- p — s- i^e + n; p -f i>e — >■ e+ -I- n. (5) 



5 



According to Bernstein [52], the rate for these reactions is: 

A = ^ (12 + 6u + u") , (6) 



where Tn is the lifetime of the neutron. In order to identify the effects of the gravita- 
tional coupling constant, we can write this lifetime as 

/i(aw,me/mp)/i . 

T~n = 5 , (7) 

nip 

where /i is a dimensionless function of the weak coupling constant, Ow, electron to 
proton mass ratio, fj,, and possibly some other mass ratios and constants, but not a 
function of Og. The freeze-out temperature is set by the equality of this rate and the 
expansion rate of the Universe: 

^ wu^-u^-6u-12 = 0, (8) 



where we have defined w through 
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and takes care of the number of relativistic species. This cubic equation can be 
solved to obtain the freeze-out temperature as a function of gravitational couphng: 

1 (ISu- + l)^/^ + (18u) + l)"'''^' + 1 
^^=0 • (9) 

After this time the most important remaining reaction is the /3-decay of neutrons, 
which continues until the formation of deuterium. Deuterium formation is delayed 
due to the large photon-to-baryon ratio. Its abundance is fixed at a temperature Tb , 
which is roughly given by the equality of the number density of photons with the 
number density of the deuterium which has been formed: 

exp(^).^l. (10) 

Here 77 is the ratio of baryons to photons and Bd is the binding energy of deuterium, 
which can be written as a dimensionless function of the strong and electromagnetic 
fine structure constants times the proton mass, Bu — /2(as, Qcm) mp c? . The function 
/2, like /i, is dimensionless and does not depend on Og. The temperature Td could 
also be converted to an age through the Friedmann equation. By the time of deu- 
terium formation neutrinos have already frozen out and electron-positron pairs have 
annihilated. Thus the ratio of the age of the Universe to the neutron lifetime becomes 



t 1 / 45 



/i/l ln(,7)2 V 32 7r3 
and the primordial helium fraction (by mass) is 



a, 



^1/2 



(11) 



>'p^2^e-^ (12) 

We have aimed at an expression for the helium abundance which is manifestly dimen- 
sionless, i.e. it only depends on dimensionless ratios of physical quantities, including 
Qg. This simple analysis leads to a primordial helium fraction of about Yp ~ 0.22, 
which is within 10 percent of the value coming from more complicated numerical BBN 
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codes. Now we can put all this together to track the effects of a possible variation of 



Og on Ypi 
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(13) 



The first term in the square brackets comes from the change in the neutron's freeze- 
out fraction and the second is due to a change in the age of the Universe. Both terms 
have the same order of magnitude (~ 10~^) and have the same sign. A higher Og 
leads to more neutrons at the freeze-out time and a lower age for the Universe, both 
of which enhance the primordial helium fraction. 

If one takes the measurement error on to be ~ 0.005, and assumes a power-law 
variation with time of the form Og oc this simple analysis shows that x should be 
less than 0.005, which is consistent with the results of other studies (see e.g. [51]). 

We end this section with a discussion of other published studies of BBN. Despite 
this literature being almost entirely dimensionful, we find that most of the papers on 
variation of G within BBN can be considered to be valid if one simply reinterprets a 
variation in G as a variation in Og. However, this is not always the case, particularly 
where extra physics is considered. Sometimes statements are made like "whenever 
you see G, interpret it as G times some quark mass m^, or a combination of G and 
Aqcd" (see e.g. [9]). But actual measurements of G have always been made using 
normal atoms, and therefore when using physical equations to derive some proposed 
time variation for G (or more properly, Og), then Og is effectively the one introduced 
in equation ([l}. In this sense, those papers which are assuming a simultaneous time 
variation for G and mp, are not valid (see e.g. [421 [9] 1. 



4 Redshift of recombination 

Cosmological recombination of hydrogen is mainly controlled by the population of the 
first excited state [53]. This is because of the high optical depth for photons coming 
from transitions direct to the ground state. The rate of recombination to this first 
excited state is given as [541 [55] 

r = „„."-9^™„e(i|)'",„(i|)^§, (14, 

Here "B" stands for baryons (i.e. protons and neutrons), "e" for electrons, q'-^-' is the 
recombination rate to the second energy level of hydrogen and -Bo is the binding energy 
of hydrogen, which is equal to a^^m^f? /2. We have assumed that all of the atoms 
are ionized. It also simplifies things considerably (without qualitatively changing the 
physics) if we ignore the mass difference of protons and neutrons, i.e. set me = mp, 
so that 

'^B = — = — Ur ■ (15) 
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This gives a rate 
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Assuming for further simplicity that Qm = 1 (this assumption could easily be relaxed 
later), the Hubble constant is 

and one can also use the equalities 



POB = , PO = PORTT— , (18) 

Hr iZR 



along with 



Putting all of these together, one can find an expression for the redshift of recombi- 
nation: 

1 I 1/2 -3 3/2 0-3/2 / ^^M \ / SIb \ /„n\ 

l + ,^o.-oc«/ a.^M^ j j . (20) 

Again we can see that an observable, which is the redshift of recombination, de- 
pends only on the dimensionless ratios of physical quantities, together with some 
other dimensionless parameters. Although this expression contains much of the essen- 
tial physics, unfortunately the numerical factor (which we neglected to write down) is 
far from correct. That is because we have not taken account of the partial hydrogen 
ionization. It would be possible to take this further by using an approximation to 
the Saha equation to correct for the ionized fraction, which would change the scalings 
with the dimensionless parameters. However, this rapidly gets complicated, and for 
accuracy one really needs the numerical solution to the relevant differential equations 
(which we will do in the next section). 

Nevertheless, the main point is that one can see an important dependence on the 
dimensionless quantities of interest. The change in the redshift of recombination is 
one of the main cosmological effects of varying Ofg. There are also other effects on 
CMB anisotropics, as we will see in the next section. 



5 Varying constants and CMB anisotropies 

So far we have written down expressions which were explicitly dimensionless, but at 
the expense of accuracy. Let us now numerically explore the effects of a varying Oem or 
ag on the CMB anisotropies. We use recfast vl.5 [56j for the recombination history 
of the Universe and CMBFAST [S7] for the calculation of CMB angular power spectra. 
It turns out that one would face some difficulties if one were to blindly dive into the 
codes and try to change "G" or add some factors of Qcm in the relevant parts. There 
are several issues which should be considered before one can promote these constants 
into dynamical parameters in the code. It turns out that this choice not only affects 
the feasibility of the problem, but also the outcome of the numerical calculation. As an 
example, CMBFAST is written so that it uses the total density contributions, {fi}, both 
for the strength of inertia in acoustic modes, and as gravitational charge. Therefore, 
if one chooses these ratios as the free parameters, in place of densities, it is much 
easier to trace the effects of a variation in the gravitational constant, because there 
is an additional factor of G wherever the code needs the Qs as seeds for gravitational 
collapse. 
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5.1 Recombination through RECFAST 

One can trace the effects of a varying fine structure constant or gravitational constant 
on the process of recombination through recfast (see e.g. |58l 1591 [60l I61p . In order 
to convert the code to a form where it can run with a different Qcm, one has to track 
all of the relevant dependencies Most of this shows up in the energy levels, 

since E oc a^^. The "Case-B" recombination rate (see Rybiki and Lightman p. 282 
[63) ) has an dependence and the complete analytic form of the 2s-ls two photon 
rate ( [641 165) ) varies as a^^- Triplet transitions of helium (discussed in [66]), case-B 
recombination and two photon transition rates for helium have the same scaling as for 
hydrogen. With this information in hand one can trace the effects of a different Ofem 
(or a time-dependent Oem) on recombination. The G dependence for recombination 
lies entirely in the Hubble constant. Putting the rates together and defining the 
dimensionless ratios, 

A=^, B = KAnp, C = KPn^, D = ^, (21) 
the Saha equation [56] will take the form 

^ ^ A - (1 - ..) ^] . (22) 

This is explicitly for hydrogen, but there is a similar form for helium too. These 
dimensionless ratios have the following dependencies on the coupling constants: 

j4 (ftom, ftg) oc Ocm ftg ''■^ ; (23) 

B (oom, Qg) OC al^ a~°'^ ; (24) 
C (acm, Qg) oc a~nj ot~°'^ ; (25) 
_D (Qcm, ag) oc ctg • (26) 

Figure[T]shows the results of a 1% increase in the fine structure constant and separately 
of the gravitational constant through the history of recombination. This is consistent 
in sign with equation pO)) - an increase in acm leads to a higher ionization fraction, 
which leads to a lower redshift for recombination. It is also noticeable that even a 
1% increase in acm can lead to more than 1% variation in the recombination history, 
while such a variation in ag does not leave any significant trace, due to the weaker 
power-law scalings above and the relatively thin last scattering surface. 



5.2 Perturbations through CMBFAST 

The effects of an alteration in the recombination history could in principle be measured 
through the CMB power spectra. This section complements section 3 in terms of 
comparing a hypothetical variation in the fundamental constants with known physics 
at a particular epoch, z ~ 1000 in this case, rather than z ^ lO'"* for BBN. As is 
discussed in [55], a different fine structure constant leaves its main imprint on the 
recombination history, and this shows up in CMBFAST through two main effects: (1) 
the derivative of the opacity, which is proportional to Thomson scattering, with an 
Oom dependence; and (2) the number density of free electrons, which is basically the 
"freeze-out" value at the end of recombination. Therefore, one can almost ignore the 
effects of a different fine structure constant after recombination. 
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Figure 1: Effects on the ionization history of 1% increases in a^m (left) and ag (right). 



The case is different for Ofg. As is discussed in [67j, a different G will lead to a 
different sound horizon at the last scattering surface and a different distance from 
this surface to us, the combination of which will change the position of the peaks of 
the power spectrum. There is also an integrated Sachs- Wolfe effect (ISW) which can 
enhance the amplitude of perturbations. This is examined in [68] and, as explained 
in [69] for the case of a constant but different G, the effects are much less important 
than they are for a time-varying G. This is partly because there is no ISW effect, but 
also because the changes on the sound horizon and distance to last scattering partially 
cancel. The results of a 1% increase in Oem and Og are shown in Fig. [2] 

In [68], the authors propose a scaling in G, (G G), which appears equivalent 

to the scaling of Qg considered here. However, we cannot confirm their results, since 
they have assumed that {To,Hq, Qb.} have the same values as in the standard model. 
This is incorrect, because these three are not independent and a change in Og will 
result in a different value for at least one of them (depending on which one we have 
chosen as our dependent variable). 

An example of time varying ag or aom results are shown in Fig. [3] We have explic- 
itly used a <x ^-o ^oo^ easy to instead imagine power-law variations in redshift, 
conformal time or some other variable. Specific ideas for time- varying constants may 
lead to specific forms for the function of time, e.g. aem(t)- Following such detailed 
predictions are beyond the scope of the present paper. However, the methodology set 
out here should be useful in testing any explicit model that is proposed. 

Our study is certainly not the first to explore the effects of variable fundamental 
constants on the CMB (see e.g. [MUZQlISS] for G and [HI [ZSl [731 [M [ZS] for aem). The 
main point here is that we have insisted on doing everything in a dimensionless way, to 
avoid finding any apparent effects which are simply a result of the choice of units. As 
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Figure 2: Effects on the CMB anisotropies of a 1% increase in (left) and (ri] 
The vertical axis is percentage. 

an example, [T^ suggests that one can search for a violation of the Strong Equivalence 
Principle using the CMB power spectra. The method, though at first sight quite 
elegant, is totally dimensional. When one starts to make it dimensionless, it becomes 
evident that the whole argument has to be reconsidered. A basic assumption is that a 
variation in G can make a difference in gravitational mass relative to the inertial mass 
of baryons. However, this is in general wrong, since one can always use the freedom of 
setting the gravitational charge and inertial mass of a given body (and only for that 
one object) equal to each other. 

Several particular ideas for exploring the variation of G within cosmology have 
been described elsewhere. Although we have not comprehensively checked all previous 
studies, we are suspicious of any in which the dimensionful quantity G is discussed 
along with other dimensionful constants, such as Ho (e.g. [771 1781 [79l 1801 181| ). 

5.3 Fisher matrix analysis 

Let us be more explicit in defining the cosmological model which describes our Uni- 
verse. In terms of the standard cosmological picture, one could define dimensionless 
parameters at some fiducial epoch. If we choose matter-radiation equality as this 
epoch, then we can define the background cosmology by giving the values of the f^s 
then, as shown in Table [T] We also need to give a quantity to set the epoch at equality, 
and we choose 9. Finally we need to provide one parameter (out of a wide range of 
possibilities) to define the time today. We have pointed out how it is possible to reach 
invalid conclusions about the variation of fundamental constants (particularly G) by 
fixing too many parameters in the cosmological model. In this section we show this in 
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Figure 3: Effects on the CMB anisotropies of time- varying constants. We have specificahy 
used Qjg (X (igft) and oc t-o oo^ (j-jght). The vertical axis is percentage. 

practice through a Fisher matrix calculation to forecast parameter uncertainties for 
future observations. 

The details of Fisher matrix analysis for CMB anisotropies are fully explained in 
[82) . Here we have used the characteristic values for the Planck satellite [H^, and 
summed over the {TT, TE, EE} channels, where T stands for temperature and E is 
the E-mode polarization. 

Specifying the value of Qa at equality and the ratio X{= flcDM/^B) suffices 
to give the correct fraction of the energy density in each sector at equality, i.e. 
{^CDM' ^B*' ^1^'''' ^7''' ^A*}' ^ fl^^ Universe. Therefore, the equality epoch is com- 
pletely defined by the set of parameters {Q']^ , X,6eq}. And the description of the 
Universe's observables will be complete by providing the time today. To (or ^o)- 

The publicly available CMB anisotropy codes are written in terms of the parame- 
ters today, rather than at a physical epoch, such as equality. To use these codes with 
variable G it is important to first ensure that the simpler situation is being considered 
properly, in which we have one value of G at recombination and another value today. 
This is done by relating the values of the parameters at equality with those defined 
today. The easiest of the set to treat is X, which is constant if neither the dark 
matter particles or the baryons can transform into anything else after equality. The 
temperature Tcq is ToQm/^h, where Q.m = S^cdm + and Q.r = fi-y -f fi^. The dark 
energy at equality is 

«A = 4- (27) 

The set of parameters chosen for the Fisher matrix calculation are: {X, 9cq, 60, t, As, n, Og, 
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Table 1: Dimensionless cosmological model parameters. We have assumed massless neutri- 
nos here, contributing to the "radiation" by the usual factor of 0.68 times the photon energy 
density. We have also assumed a flat background, and calculated values and uncertainties 
using the Markov chains from the WMAP 6-parametcr fits [37] . 

Here, r is the optical depth of the Universe after reionization, As is the amplitude 
of the power spectrum of scalar perturbations (which could easily be made dimen- 
sionless) and n is its power-law slope. One could also perform the same calcula- 
tion in a simplistic and incorrect way by choosing the wrong set of parameters, 
{/lo, ficDM, f^B, T, ^B, n, Q!g, Oc}, that are both time dependent and inter-dependent 
on each other, and by over-constraining the observation epoch, i.e. ignoring Equa- 
tion [2l In order to perform a fully consistent analysis, one should also take care to 
normalize the power spectra to the initial conditions and not to the large angle CMB 
anisotropics today; CAMB[84] was chosen for the Fisher matrix analysis, since one 
can choose either option for normalization. 

The 1(7 error bars for the two different methods are contrasted in Table [2] and 
the correlations among the parameters are shown in Fig. [l] The conclusion is that 
the error bar on Qg changes by approximately a factor of 2 between the two different 
approaches, and Qg shows a slightly stronger correlation with other parameters in the 
correct approach compared to the incorrect one. This underscores the need to be 
careful when considering the effect of variable constants on cosmological observables. 

6 Other cosmological observables 

We have already considered the effects of a variable Qg or Qcm on BBN and CMB 
anisotropics, with the main emphasis on making things dimensionless. One can extend 
this approach to other observables in cosmology, such as weak lensing, baryon acoustic 
oscillations (BAO), large-scale structure, ISW effect and the use of supernovae as 
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Parameter la values 

correct method incorrect method 

hi(acm) 0.0058 0.0056 

hi(ag) 0.031 0.062 

ln{9eq) 0.0096 0.04* 

ln(6lo) 0.0016 0.025* 

X 0.10 0.17* 

^ X 10^ 0.11 0.105 

n 0.0074 0.0073 

T 0.021 0.021 

ho 1.67* 1.38 

Qb 0.0017* 0.00091 

f^cDM 0.0056* 0.0032 



Table 2: Comparison of the la error bars of the cosmological parameters of the two different 
methods described in Section r5.3l for the Fisher matrix analysis. The first approach involves 
the cosmological model in an explicitly dimensionless way, while in the second approach 
G is simply added as a parameter to the usual model. A star on a value shows that the 
parameter was not among the set of parameters for Fisher matrix calculation in the relevant 
method. 

standard candles. We will avoid discussing the use of high redshift quasar spectroscopy 
to constrain Qcm and fi, since that topic is extensively covered elsewhere (see references 

in m- 

6.1 Baryon acoustic oscillations 

BAO follows the same basic physics as CMB anisotropics (see e.g. [85]), so there 
should in principle be the same kind of Og dependence. The important timescales 
for BAO are the equality epoch, Zeq and the drag epoch, Zd- The first, Zcq, is purely 
determined by the initial conditions of the Universe, and the time today, setting the 
epoch when = i^R- The drag epoch is defined as the time when baryons are freed 
from the Compton scattering of photons, and therefore Zd has the same Og or aem 
dependence as Zr. These special redshifts set the scaling conditions for BAO. After 
the drag epoch, one should solve for the matter transfer function, where this function 
is basically the same as the CMB transfer function in terms of Og or Qem dependence. 
It seems clear therefore, that if one wanted to use BAO to constrain the variation of 
Og, it would be straightforward to follow the same procedure we discussed in section 5. 

6.2 Gravitational Lensing 

In gravitational lensing (see [86] for an earlier study) there is a danger of becoming 
confused by dimensions, since both the estimate of the lensing mass and the curvature 
depend on G, and they are mixed together in Qg. The simple case of an Einstein ring 
is illuminating. Here the lens and the source are colinear and in Euclidean space 
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In(ag) A n i Ha^) A n t 



Figure 4: Comparison of the eorrelation among selected parameters for the two different 
methods described in Section 15.31 for the Fisher matrix analysis. The diagram on the left 
corresponds to the correct method, while the panels on the right show the results of the 
naive (and incorrect) approach. 



the radius of the Einstein ring is the geometric mean of the Schwarzschild radius of 
the lens and the distance to the lens. This means that for a relatively nearby lens, 
at distance d, the lensing angle is Oe ~ \/ GM/c^d, and if we estimate the mass 
through measuring a velocity dispersion over a radius r, then 0e ~ /c^ , 

where = r/d is the apparent angular size of the object. Viewed this way we see 
that it is hard to use lensing to measure the strength of gravity, because the obvious 
dimensionless observables leave no G dependence! 

Let us look at this in a little more detail. In cosmology one can work out the 
lensing angle to be (e.g. [ST]) 



where dLS is the distance from the source to the lens, and ds and cJl are, respectively, 
the distances from us to the source and to the lens. However, this seemingly innocent 
dimensionless equation is not useful for any experiment designed to measure a change 
in the gravitational constant; this equation, or any equation for lensing, should be 
turned into an equation which only has Qg or Qcm dependence before it can be used 
as a test of fundamental constant variation. One can work out all of the distances in 
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the above equation and end up with the following relation: 



with A 



Be 




1 + 



1 + ZL 



(29) 



(30) 



Here and zs are, respectively, the redshift of the lens and the source and A^l and 
A^s are defined as: 



where x can be either S or L. 

Speaking in more general terms (and again thinking about communicating with 
another civilization), we can see two possibilities. The first is that one could imagine 
choosing a lens which consists of a fixed number of particles, so that one knows the 
value of the pure number M/nip. That may be interesting from a philosophical point 
of view, but is hardly related to how we carry out observations in cosmology. A second 
idea is that perhaps one could perform a statistical survey over many galaxy lenses, 
measuring statistics which depend on a characteristic galaxy mass Mgai. If that mass 
depends on fundamental constants through a cooling argument (e.g. [88]), then it 
may be that in a lensing survey the observables depend on the number Q^j^Qg ^ ^ ^^"^ . 
Since this would be different in universes with different values of the constants, then 
this is potentially measurable, and hence it might be possible to constrain a redshift 
dependence of lensing observables. 

These ideas do not seem particularly practical, and so we are left unclear about 
whether gravitational lensing could ever be used to constrain the time variation of 
ag. Things will presumably be unambiguous in explicit self-consistent models which 
contain a variable strength of gravity, and we leave this topic for future studies of 
specific models. 

6.3 Large-scale structure and supernovae 

There are many other astrophysical phenomena which can in principle be used to test 
the variation of fundamental constants (see e.g. |13l \35\ ). Large-scale structure and 
various measures of the power spectrum of galaxy and matter clustering [89] , will also 
have dependence on a^. It will again be important to ensure that this dependence is 
dimensionless, and that cosmological parameters are not over-constrained when doing 
so. In other words, one needs to realise that changing Og can effectively change the 
epoch at which we live. 

Supernovae have also been very useful as approximately standard candles in cos- 
mology, and such studies can also be adapted to constrain a combination of fundamen- 
tal constant variations [90U91ll92]|93| . Like other stellar sources, there is a dependence 
on Og^^'^ in the evolutionary timescale (see e.g. |94l 1951 [96l 1971 198| ) . and hence the use 
of supernovae, combining the standard candle property with luminosity distance, will 
involve a different combination of fundamental constant variations than for BAO and 
lensing studies. A combination of cosmological probes will therefore enable variations 
among different parameters to be distinguished. 

Finally, we note that among the different arenas within astrophysics, cosmological 
studies have the potential of probing not only the time dependency of ag, but also 




1 + Zx 



(31) 
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any scale dependence. It is important to work out an appropriately consistent theory 
to describe these phenomena within the context of a spatially variable fundamental 
constant. Although we have not considered such models here, it will surely still be 
important to avoid dimensional quantities. 

7 Conclusions 

We are not advocating that all of cosmology henceforth should be represented in 
dimensionless forms. However, we would say that care has to be taken when dealing 
with variation of physical constants in cosmology. This is particularly because of 
the cosmic time ambiguity, which could lead to oiier-constraining the cosmological 
parameters, especially when G is involved. 

Definite physical mechanisms for the variation of fundamental constants will lead 
to specific forms for ag(f) etc. We believe that the proper place to start investigating 
such theories is to make sure that there is a robust basis for comparing cosmologies 
with different constant values of the constants. Only then can one effectively deal 
with time-variable quantities. 

A time-dependent (or space-dependent) Qg is not consistent with General Rela- 
tivity. One explicit framework for accommodating such a variation are the scalar- 
tensor theories of gravity. Some studies have already investigated CMB anisotropies 
and other cosmological constraints in the simplest form of scalar-tensor model (e.g. 
[99lllOOII101] l. the so-called Brans-Dicke theory |102] . While this version of a scalar- 
tensor theory seems to be already ruled out by solar system experiments [103| . there 
are more general scalar-tensor theories which pass the solar system tests |104| . and 
might be promising avenues of exploration (e.g. [T] ?]). These may provide analogues 
for investigating the variation in physical constants proposed in some brane-world 
scenarios. In studying the empirical tests of such models we recommend keeping ev- 
erything dimensionless in order not to be misled by apparent variations that may be 
unmeasurable. 
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